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The impact of recent multi-loop calculations on precise determinations of charm- and bottom-quark masses and the 
strong coupling constant is discussed. 



QCD has developed from a qualitative model for strong interactions into a quantitative theory with precise pre- 
dictions for a multitude of observables. This development, which has taken place in particular during the past 
two decades, is due to significant improvements in our theoretical understanding, in calculational techniques and 
in improved experimental results, ranging from low energy studies in r-lepton decays through precise cross section 
measurements in electron-positron annihilation up to measurements of the Z-boson decay rate into hadrons with 
permillc accuracy. Many of these inclusive observables are related to essentially the same object in quantum field 
theory, namely the absorptive part of the current-current correlator. Restricting the discussion to measurements at 
high energies or to properly chosen integrals, these observables can be evaluated in perturbation theory and lead to 
precise determinations of the fundamental parameters of QCD, quark masses and the strong coupling constant. 

1 . Quark masses 

A detailed analysis of TOc a-nd mt, based on the ITEP sum rules [1] has been performed several years ago Q. During 
the past years new and more precise data for a{e~^e~' —>■ hadrons) have become available in the low energy region, in 
particular for the parameters of the charmonium and bottomonium resonances. Furthermore, the error in the strong 
coupling constant as{Mz) which enters this analysis has been reduced. Last not least, the vacuum polarization 
induced by massive quarks has been computed in four-loop approximation; more precisely: its first derivative at 

— 0, which corresponds to the lowest moment of the familiar i?-ratio has been evaluated in [1, 0|- Based on these 
developments a new determination of the quark masses has been performed in Ref. 5] . More recently also the second 
moment has been calculated fo^. 

The extraction of toq from low moments of the cross section (T(e+e^ QQ) exploits its sharp rise close to the 
threshold for open charm and bottom production. By evaluating the moments 



with low values of n, the long distance contributions are averaged out and A4n involves short distance physics only, 
with a characteristic scale of order -Ethioshoid = 2toq . Through dispersion relations the moments are directly related 
to derivatives of the vacuum polarization function at — 0, 



which can be evaluated in perturbative QCD (mg = rriQ^fi) is the MS mass at the scale /z). The perturbative series 
for the coefficients Cn in order was originally evaluated up to rt = 8 in Ref. J6[ , and to "arbitrary" high order 
in 0, S|- The four- loop contributions to Co and Ci were evaluated in Refs. [j], those to C2 in Q. Combining 
Eqs. ([1]) and ([2]) the quark mass can be extracted. At this point it should be emphasized that the relative weight 
of resonances and continuum is quite different in the experimental moments. Furthermore, low moments are less 
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Figure 1: Comparison of recent determinations of mc(3GeV) and mt{mb) 



sensitive to non-perturbative contributions from condensates, to the Coulombic higher order effects, the variation of 
/X and the parametric dependence. For n = 1: 



mc{3 GeV) - 0.986(13) GeV 



(3) 



The moment with n = 2 is less sensitive to data for R{s) from the continuum region above 5 GeV, where experimental 
results are scarce and the aforementioned theory uncertainties are still relatively small. The agreement between 
n = 1 and n = 2 (mc(3GeV) = 0.976(16) GeV), together with the nice convergence with increasing order in as can 
be considered as additional confirmation of this approach. 

Instead of measuring the moments A^^'^p in e~^e~ annihilation they can also be determined in lattice simulations. 
This approach has recently been pioneered in using the Highly Improved Staggered Quarks (HISQ) discretization 
of the quark action in combination with four-loop perturbative results BQSHllll. The final result, mc(3 GeV) 
0.986(10) GeV corresponds to a scale-invariant mass mc(mc) = 1.268(9) GeV and is in excellent agreement with the 
determinations based on e+e^ data. 

The approach based on e^e~ data is also applicable to the determination of to;,. The three results based on n = 1, 2 
and 3 are of comparable precision. The relative size of the contributions from the threshold and the continuum region 
decreases for the moments n — 2 and 3. On the other hand, the theory uncertainty is still small. Therefore the 
result from n = 2 was taken as the final answer despite the fact that C2 was not yet known. The result, 
™fc(10 GeV) = 3.609(25) GeV, corresponds to mhipib) = 4.164(25) GeV. The recent evaluation Q of C2 has lead to 
a decrease of the central value by 2 MeV and a reduction of the error from 25 MeV to 19 MeV. A comparison of a 
few selected m;,- and mc-determinations is shown in Fig. [TJ 



2. The strong coupling constant 

One of the most precise and theoretically safe determination of as is based on measurements of the cross section 
for electron-positron annihilation into hadrons [l^ . These have been performed in the low-energy region between 
2 GeV and 10 GeV and, in particular, at and around the Z resonance at 91.2 GeV. Conceptually closely related is 
the measurement of the semileptonic decay rate of the r-lepton, leading to a determination of as at a scale below 
2 GeV 13[. The perturbative expansion for the ratio R{s) = a{e^e~ — > hadrons)/(T(e+e^ /i+/i^) in numerical 



form is given by 

R = 1 + as + (1-9857 - 0.1152n/)a2 + (-6.63694- 1.2001371/ - 0.00518n^) a 

V 



(-156.61 + 18.77n/- 0.797471^ + 0.0215^)0^ . (4) 



Here as = as/ir and the normalization scale jj? — s. The corrections are conveniently classified according to their 
power oi Uf, with Uf denoting the number of light quarks. The a'^n^j^ term is part of the "renormalon chain", the 
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next term of order a'^n'j was evaluated in the complete five- loop calculation has been performed in 15 1. 

Let us now move to the analysis of present data for e+e^ annihilation and r decays. Measurements of R{s) at lower 
energies, with their correspondingly larger values of as, are in principle more sensitive to as{Mz) if the same relative 
precision could be obtained. At present, however, the systematic experimental error of 2% is a limiting element 
for a competitive measurement. The final result of a recent analysis [l6| ai^-'(9GeV) — 0.182lQ g25 represents the 
combined information on the strong coupling from R measurements in the region between 3 GeV and the bottom 
threshold. In contrast to the situation below 10 GeV the extraction of as from Z-decays is affected by the terms. 
The 0{al) analysis of the electroweak working group [l^l is shifted 15| by 6as{Mz) = 0.0005 and one finds 



as{Mz)^^^^ = 0.1190 ± 0.0026°'^P . (5) 

The theory error can now safely be neglected. 

Higher orders are of particular relevance in the low-energy region, for example in r decays. The correction from 
perturbative QCD to the ratio T(t — > hadronss^o + '^r)/r(T I + i?i + v^) is given by the factor 

'"'^ ds s V f 2s 







which can be evaluated in Fixed Order perturbation theory or with "Contour Improvement" as proposed in 1^ 19 1: 

5^^ = fls + 5.202 al + 26.366 + 127.079 , 5^^ = 1.364 + 2.54 al -f 9.71 4- 64.29 . (6) 

Starting from ^q'^p = 0.1998 ± 0.0043oxp one obtains 151 a™(Ai^) = 0.322 ± 0.004exp ± 0.02 and a^^{Mr) = 
0.342 ± 0.005exp ± 0.01, respectively. The second uncertainty corresponds to a change in the renormalization scale 
/Lt^/M^ between 0.4 and 2. 

As stated above the theory error for ag from Z decays is small compared to the experimental uncertainties. For r 



decays the difference between FORT and CIPT must be considered as irreducible uncertainty [15| : 

as {Mr) = 0.332 ± 0.005oxp ± 0.015th ■ (7) 
Applying four- loop running and matching 21, 22, 2^, 24 1 (with negligible error from the evolution from Mr to Mz)'- 



as(Mz) = 0.1202 ± 0.0006oxp ± 0.0018th = 0.1202 ± 0.0019 . (8) 

The shifts in as{Mz) from Z- and r-decays, are opposite in sign and move the values in the proper direction, 
decreasing, thus, the current slight mismatch between two independent determinations of as- The two results are in 
remarkable agreement and can be combined to one of the most precise and presently only N^LO result: 

as{Mz) = 0.1198 ± 0.0015 . (9) 

As discussed above, as from r-decays is strongly affected by theory uncertainties. This is reflected in three recent 
publications which all are based on the same set of t data and the new five-_loop results, which, however, arrive at 
significantly different estimates of higher order terms and hence of as. In 25] it is argued that FOPT exhibits a 
poorly "convergent" series, in contrast to CIPT, where a subset of higher order contributions (so-called 7r^-terms) is 
automatically summed. This analysis leads to as (A/,-) — 0.344 ± 0.005ea;p ± 0.007fh. The small theory uncertainty 
is a consequence of the (artificially) reduced /i-dependence and the restrictive assumptions about higher orders in 



perturbation theory. The opposite viewpoint has been advocated in 26[ where an explicit and plausible behaviour 
of the perturbative series is modelled. Strong cancellations are observed between the aforementioned 7r^-terms and 
those terms which can only be obtained from a concrete higher order calculation. 

All these papers are based on similar assumptions about non-perturbative power-suppressed contributions. These 



assumptions have been questioned in 27[, employing different weight functions in order to suppress these poorly 
determined terms with dimension D > 8. A significantly smaller result, asiMr) = 0.3209 ± 0.0046 ± 0.0018, is 
obtained within CIPT. In total the spread among the different results is covered reasonably well by eq. ([7]). 
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Last not least it is instructive to compare eqs. (O and ([8]) with as determinations based on lattice simulations. 
The most recent results from a precision simultaion by the HPQCD collaboration [2^ as{Mz) — 0.1183 ± 0.0007 
and from a study [2^ based on essentially the same lattice data and an alternative perturbative analysis, as{Mz) = 
0.1192 ± 0.0011 are in mutual agreement, as should be expected because they are based on a similar set of data. 

Remarkably enough, a completely different approach [lO|, which exploits the lattice simulation of the pseudoscalar 
correlator and which has been mentioned already above leads to as{Mz) = 0.1174 ± 0.0012 again well consistent 
with the other determinations. 
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